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CHAPTER  I 


INTRODUCTION 


Various  methods  have  been  attempted  to  predict  the  cold  circuit 
characteristics  of  coupled  cavity  circuits  for  traveling  wave  tubes.  These 
methods  include  approximate  field  analyses  [l]  and  lumped  element  equivalent 
circuits  [23,  None  of  these  methods  can  predict  the  circuit  characteristics 
to  a  high  degree  of  accuracy.  The  result  is  that  the  circuit  design  procedure 
includes  physically  building  and  measuring  a  "first-cut"  circuit,  and  making 
modifications  to  obtain  the  final  product.  These  modifications  include  both 
engineering  and  machining  costs. 

With  the  advent  of  large  digital  computers  with  low  price/performance 
ratios,  and  the  development  of  numerical  techniques  for  efficiently  handling 
the  massive  amounts  of  data  storage  necessary  for  self-consi stent ly  solving 
Maxwell’s  equations  with  sources.  It  has  become  financially  feasible  to  make 
use  of  computer  codes  which  perform  three-dimensional  "numerical  experiments." 
Codes  of  this  type  have  been  used  for  a  number  of  years  in  Plasma  Physics, 
and  particularly  for  EMP  calculations. 

This  study  makes  use  of  a  three-dimensional  plasma  simulation  code 
called  SOS  (Self -Optimized  Sector)  [3],  developed  by  Mission  Research  Corporati 
By  making  certain  modifications  to  the  source  models  and  using  the  output 
data  of  the  code  as  input  to  a  Fast  Fourier  Transform  routine,  it  was  possible 
to  perform  numerically  the  Resonance  tests  by  which  a  circuit's  Brillouin 
plot  is  typically  found. 


Chapter  II  will  introduce  the  code  and  its  method,  which  is  a  great 
contrast  to  the  eigenvalue  approach  used  in  codes  such  as  Superfish  [4]. 
Chapter  III  discusses  the  various  source  methods  which  were  used.  In 
Chapter  IV  various  test  cases  are  discussed  which  establish  the  validity  of 
the  method.  Chapter  V  gives  the  results  of  a  staggered-slot,  coupled  cavity 
circuit,  and  compares  these  to  measured  values.  The  problem  of  resolution 
and  an  apparent  numerical  instability  are  discussed  In  Chapter  VI,  with  an 
eye  towards  cost  effectiveness.  Chapter  VII  gives  a  summary  of  results  and 
recommendations  for  further  study.  In  Appendix  A,  the  second  order  accuracy 
of  the  non-uniform  grid,  centered  difference  derivative  method  is  proved. 
Appendix  B  gives  a  derivation  of  the  Discrete  Fourier  Transform  pair.  The 
equations  necessary  for  fitting  a  cubic  least  squares  with  boundary  condi¬ 
tions  (zero  first  derivative)  are  developed  in  Appendix  C. 


In  component  form,  the  two  curl  equations  become 
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Where  3j  denotes  the  partial  derivative  with  respect  to  the  i*h 
coordinate,  and  the  coefficients,  hj,  are  the  differential  length  elements 
in  the  coordinate  system  of  interest,  e.g. 

in  rectangular  coordinates:  hj=1;  h2r1;  hs*l 

in  cylindrical:  hi = 1 ;  h2=r;  h3  =  1 

and  in  spherical:  hj=l;  h2=r,  h3=rsin0. 

One  can  apply  various  finite  difference  schemes  to  these  equations 
based  on  accuracy,  convergence,  and  cost  considerations.  505  contains 
four  finite  differencing  algorithms  from  which  the  most  appropriate  for 
a  given  application  may  be  chosen.  Of  these,  two  deal  with  regions  of 
temporally  and  spatially  varying  conductivity  and  need  not  be  considered 
for  our  application.  A  third  algorithm  uses  the  centered  difference  formula 
for  the  spatial  derivatives,  but  is  implicit  in  the  time  domain  of  the 
magnetic  field.  That  is,  a  "guess"  at  the  future  value  of  the  magnetic 
field  is  used  to  calculate  the  current  electric  field.  This  electric 
field  is  then  used  to  calculate  the  next  magnetic  field  value  which  is 
compared  to  the  guess.  By  adjusting  the  guess  appropr lately,  the  algorithm 
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converges  to  the  desired  field  values.  The  fourth  method  is  a  straightforward 
centered  difference.  In  this  approach,  the  electric  field  is  found  exclusively 
from  the  previously  calculated  magnetic  field  values.  It  is  therefore  explicit 
so  that  no  iterative  procedure  is  needed. 

Although  the  third  approach  will  yield  more  accurate  results,  it  is 
clearly  a  more  expensive  method.  For  this  reason,  the  explicit,  centered 
difference  algorithm  was  used  for  calculations  in  this  report.  This  method 
can  be  shown  (see  Appendix  A)  to  yield  second  order  accuracy  to  the  true 
derivative,  while  retaining  a  fairly  flexible  grid  spacing  requirement. 

SOS  contains  the  option  to  sample  any  of  the  field  components  at  a 
fixed  point  in  space  as  a  function  of  time.  This  is  accomplished  by  reading 
the  field  amplitude  at  discrete  points  in  time  of  fixed  interval,  and  writing 
the  result  to  an  output  file.  It  is  straightforward  then,  to  apply  a  Fourier 
Transform  to  this  data  so  as  to  obtain  the  frequency  characteristics  of  the 
device  under  test.  To  this  end,  we  introduce  the  Discrete  Fourier  Transform 
pair  [5] 
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where  N  is  the  number  of  samples  taken.  (See  Appendix  B  for  a  derivation 
of  the  DFT  pair.) 


If  we  define  tc  as  ttie  sampling  Interval,  then  we  have 


t  = 


mtc 

TT 


and 


f  = 


Nt 


Notice  that  Nt0  is  the  total  time  through  which  the  fields  have  been  tracked 
in  our  numerical  experiment. 

The  cost  (CPU  time)  of  a  simulation  is  determined  primarily  by  two 
factors:  the  number  of  grid  points,  and  the  number  of  time  steps.  Obviously, 
we  want  to  minimize  both  in  order  to  keep  the  cost  down.  The  first  is  limited 
by  the  requirement  to  fully  define  the  physical  details  (e.g.  corners)  of  the 
structure.  The  number  of  time  steps  is  more  complicated,  as  it  is  a  function 
of  several  factors. 

The  spacing  in  frequency  space  is  determined  by  the  total  time  of 
the  run,  i .e. 

f. 


1 

Nt0 


Therefore,  for  a  desired  value  of  f0,  we  can  keep  N  small  by  making  t0  large. 
However,  in  order  to  have  stability,  tQ  is  bounded  by  the  Courant  Stability 
Criterion  [6] 


where  the  hj  are  the  coefficients  of  the  differential  length  elements  defined 
above,  and  the  6xj  are  the  differential  length  elements.  For  example,  in 
cylindrical  coordinates 
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The  final  factor  to  be  considered  is  that  in  order  to  use  a  Fast 
Fourier  Transform,  fQ  and  tQ  must  be  chosen  so  that  N  is  an  integer  power  of 
two.  This  last  consideration  did  not  pose  any  problem  in  this  study  because 
if  the  "power  of  two"  restriction  was  inconvenient,  the  FFT  could  be  replaced 
by  a  Discrete  Fourier  Transform  (not  fast).  The  extra  cost  incurred  by  this 
change  was  insignificant  in  comparison  to  the  cost  of  running  the  plasma 


CHAPTER  III 
SOURCE  CONFIGURATION 


The  sources  available  in  the  code  (SOS)  were  not  suitable  for  our 
purposes,  so  the  approach  used  was  to  manipulate  one  component  of  the 
electric  field  at  various  points  in  space.  This  was  accomplished  by 
adding  the  desired  source  (e.g.  a  point  source  oscillating  at  a  desired 
frequency,  or  a  point  impulse  function)  to  the  existing  field  value  as  a 
"correction"  term. 

This  approach  allowed  for  a  good  deal  of  freedom  In  coupling  to  the 
desired  resonance.  In  fact,  it  was  possible  to  change  the  coupling  separately 
in  physical  and  frequency  space.  By  using  an  impulse  as  the  source,  all  fre¬ 
quencies  were  excited.  If  the  impulse  was  tailored  in  space  to  fit  the 
desired  field  pattern  (which  is  known  in  the  case  of  a  simple  cavity)  then 
the  code  no  longer  sees  this  as  an  impulse,  but  rather  as  initial  conditions. 
This  is  because  we  are  using  the  algorithm  (see  Chapter  II)  which  is  explicit, 
and  our  source  is  being  artificially  added  to  the  existing  fields.  When  we 
adjust  the  field  values  at  all  points  in  the  volume  simultaneously  to 
exactly  match  one  mode,  then  there  are  no  steps  In  the  field  in  physical 
space,  and  the  code  does  not  see  the  step  in  time,  but  rather  treats  It 
as  if  that  mode  had  always  been  there  and  we  get  a  single  spike  in  the 
frequency  spectrum. 

Regardless  of  what  type  of  source  is  used,  the  energy  going  into 
the  system  will  try  to  resonate  in  the  natural  modes  of  the  system.  Most 


of  the  energy,  in  fact,  should  go  into  the  natural  mode  which  we  are  the 
closest  to,  with  respect  to  frequency  and  field  shape. 

Since  the  field  shape  is  known  exactly  in  a  simple  cavity,  and 
approximately  in  a  reentrant  cavity,  initial  conditions  proved  to  be  a  good 
way  to  find  the  primary  resonance.  Having  a  source  point  oscillating  at  a 
frequency  close  to  the  desired  resonance  also  provided  good  results. 

Figures  3.1  and  3.2  show  the  frequency  spectrum  of  a  simple  and 
reentrant  cavity,  respect i vely,  in  which  a  point  source  oscillating  at  a 
frequency  close  to  the  primary  resonance  was  placed.  Notice  that  the 
shape  of  the  curve  in  each  case  is  basically  that  of  a  peak  at  the  primary 
resonance,  with  the  source  and  higher  order  resonances  superimposed  upon  it. 
Figure  3.3  shows  the  spectrum,  after  a  beam  hole  was  introduced.  Notice 
that  now  the  shape  of  the  curve  has  changed  significantly  and  there  is 
energy  at  zero  frequency  which  is  comparable  to  the  amount  at  the  primary 
resonance.  After  introducing  the  slot,  while  still  using  a  point  source, 
Figure  3.4  shows  the  zero  frequency  amplitude  to  be  considerably  above 
the  primary  resonance.  In  this  latter  case,  if  we  look  at  higher  frequencies 
as  in  Figure  3.5,  we  see  that  the  code  is  predicting  an  incredible  amount  of 
energy  going  into  some  higher  order  mode  (about  46  GHz).  This  high  frequency 
problem  is  not  true  for  the  "no-slot"  case,  as  shown  in  Figure  3.6. 

From  these  figures,  it  is  evident  that  the  predicted  spectrum  is  not 
physically  realistic  for  either  the  case  of  a  half-cavity  with  a  beam  hole, 
or  a  half-cavity  with  a  beam  hole  and  80°  slot.  Three  possible  methods  of 
improving  the  simulation  come  to  mind  immediately.  First,  the  number  of 
grid  points  might  be  increased  to  improve  the  accuracy  of  numerical 
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derivatives.  Second,  the  source  might  be  placed  at  the  resonant  frequency 


(I.e.  Improve  the  coupling  In  frequency  space).  Finally,  the  source  might 


be  tailored  to  fit  the  field  shape  of  the  desired  mode  (i.e.  Improve  the 


coupling  in  physical  space). 


The  first  suggestion,  improving  the  grid,  did  change  the  results. 


but  did  not  solve  the  problems  mentioned  above.  This  can  be  seen  by 


comparing  Figure  3.7  to  Figure  3.5.  The  second  suggestion  amounts  to  knowing 


the  solution  to  find  the  solution,  so  it  is  not  usable.  The  third  idea 


is  useful,  because  we  know  the  desired  mode  has  a  concentration  of  electric 


field  across  the  ferrules,  so  it  seemed  reasonable  to  change  from  a  single 


point  source  to  a  ring  of  sources  located  at  all  grid  points  between  the 


ferrules.  This  method  solved  the  anomaly  of  having  large  amplitude  at 


zero  frequency,  as  shown  in  Figures  3.8  and  3.9.  However,  it  did  not 


solve  the  high  frequency  problem  for  the  cavity  with  a  slot,  as  shown  in 


Figure  3.9. 


With  these  improved  results,  it  seemed  good  to  observe  the  effect 


of  our  other  two  suggestions.  Figures  3.10  and  3.11,  with  a  finer  grid, 


may  be  compared  to  Figures  3.8  and  3.9.  Figures  3.12  and  3.15,  in  which  the 


source  is  placed  at  the  resonant  frequency,  may  be  compared  with  Figures 


3.8  and  3.9.  In  each  case,  the  change  does  not  affect  the  result  of 


interest,  which  is  the  detection  of  the  primary  resonance  and  its  location. 


Lwy 

m 

m. 


m 

4 

m 


as 


3R0I 


f 

w 

Ip 


— j.204 


TM010  RESONANCE 
(4.52  GHz) 


1.0*  SIMPLE  CAVITY 


6  7  8  9  10  11  12  13  14 

F  (GHz) 


FIGURE  3.1  -  Amplitude  vs.  Frequency  in  a  simulation  of  a  simple 
cavity  in  which  the  source  was  a  point,  oscillating  just  below 
(in  frequency)  the  primary  resonance.  Note  that  the  shape  of 
tho  curve  is  primarily  a  peak  at  the  desired  resonance,  with 
the  source  and  higher  order  resonances  superimposed  uDon  it. 
Limensions  are  inches. 
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FIGURE  3.2  -  Amplitude  vs.  Frequency  In  a  simulation  of  a 
reentrant  cavity  of  indicated  dimensions.  Again,  a  point 
source  was  placed  just  below  the  desired  resonance,  whi-ch 
was  known  from  Moreno's  curves  [8] .  The  curve  is  not  as 
smooth  as  in  the  previous  figure,  but  the  primary  resonance 
is  predicted  by  the  code  to  be  just  where  we  expect  It. 
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FIGURE  3.3  -  A  reentrant  cavity  after  the  inclusion  of  a  beam 
hole.  The  extension  of  the  beam  hole  is  to  account  for  the  wall 
thickness  between  cavities.  In  this  case,  no  predicted  value  for 
the  resonance  was  available,  but  as  can  be  seen  from  the  figure, 
the  source  was  placed  slightly  above  the  resonant  frequency. 


CAVITY  WITH  BEAM  HOLE  AND  80*  SLOT 


SECTION  A-A 
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FIGURE  3.4  -  The  final  step  to  modeling  a  half  cavity  is  the 
inclusion  of  the  coupling  slot.  This  lowered  the  resonant 

frequency  slightly,  and  changed  the  overall  shape  of  the  curve 
dramatically.  It  should  be  noted  that  this  step  makes  the 
problem  a  truly  three-dimensional  one. 
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FIGURE  3.5  -  This  figure  is  the  same  as  Figure  3.4,  but  it  shows 
all  of  the  data  obtained  from  the  FFT  (not  just  0-14  GHz  as  in 
the  previous  figure).  Two  properties  of  this  curve  are  somewhat 
disturbing.  First  the  large  amplitude  at  zero  frequency,  and 
second  the  large  amplitude  at  high  frequency. 


FIGURE  3.7  -  This  is  the  same  structure  that  was  simulated  in 
Figure  3.5,  but  with  a  better  grid  (i.e.  more  grid  points). 
Improving  the  grid  has  clearly  not  solved  either  of  our  problems, 
in  fact  it  has  made  the  zero  frequency  amplitude  considerably 
I arqer . 


CAVITY  -  NO  SLOT 
RING  ACROSS  FERRULES  AT 
3.0  GHz 
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FIGURE  3.8  -  By  changing  the  source  from  a  single  point  to  a 
ring  of  points  in  the  0-Z  plane  corresponding  to  the  ferrule 
(i.e.  improving  the  coupling  to  the  desired  mode  in  physical 
space),  the  zero  freauency  problem  of  Figure  3.6  was  eliminated 
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CAVITY  -  80°  SLOT 
RING  ACROSS  FERRULES  AT 
3.0  GHz 
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FIGURE  3.9  -  With  the  inclusion  of  the  improved  source,  the 
zero  frequency  amplitude  was  reduced  to  a  reasonable  value  for 
the  cavity  with  the  slot.  It  is  evident,  however,  that  this 
did  not  help  the  problem  at  high  frequencies. 


iURE  3.10  -  By  comparing  this  to  Figure  3.8,  the  results  of 
ng  a  finer  grid  with  the  improved  physical  space  coupling  may 
seen  for  the  cavity  with  no  slot. 
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CAVITY  -  80*  SLOT 
RING  ACROSS  FERRULES 
IMPROVED  GRID 


F  (GHz) 

FIGURE  3.11  -  This  result  is  a  sort  of  hybrid  of  Figures  3.7  and 
3.9.  The  ring  source  at  the  ferrule  has  kept  the  zero  frequency 
amplitude  low,  as  in  Figure  3.7,  but  the  improved  grid  has  intro 
duced  some  other  resonances  which  were  evident  in  Figure  3.9. 


CAVITY  -  NO  SLOT 
RING  ACROSS  FERRULES  AT 
4.0  GHz 


FIGURE  3.12  -  Placing  the  source  on  the  frequency  found  in 
Figure  3.8  to  be  the  resonant  value,  the  curve  is  very  clean 
as  one  would  expect  since  we  are  coupling  almost  perfectly  to 
a  single  mode. 
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CAVITY  -  80°  SLOT 
RING  ACROSS  FERRULES  AT 
3.98  GHz 
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FIGURE  3.13  -  Placing  the  source  on  the  resonance  does  not  have 
the  same  effect  when  our  simulation  includes  a  slot.  It  seems 
clear  that  the  high  frequency  problems  are  a  numerical  anomaly 
which  are  not  physically  realistic  and  should  not  concern  us  as 
long  as  they  do  not  affect  our  ability  to  accurately  determine 
the  true  resonance  of  our  structure. 
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CHAPTER  IV 
RESULTS 


|  As  illustrated  in  the  previous  chapter,  the  approach  used  in  this 

’ 

I 

j  project  was  to  check  the  accuracy  of  the  code  (SOS)  and  the  FFT  routine 

!  against  simple  test  cases,  and  work  towards  the  desired  structure.  The 

obvious  choice  for  a  first  test  case  was  a  right-circular-cylindrical 
cavity.  There  are  two  reasons  why  this  was  a  good  choice.  First,  the 
results  could  easily  be  compared  to  the  theoretical  solution  [7],  and  second, 
a  relatively  small  number  of  grid  points  accurately  represented  the  structure 
so  that  test  runs  on  the  code  were  inexpensive. 

Accurate  results  were  obtained  for  the  simple  cavity  case  (less  than 
1$  error),  so  the  next  logical  structure  was  a  reentrant  cavity.  A 
reentrant  cavity,  which  had  a  variable  inner  conductor  height,  was  available 
lor  measurement.  The  calculated  (by  505)  and  measured  values  are  shown 
versus  inner  conductor  height  in  Figure  4.1.  In  all  cases,  the  measured 
value  agreed  with  the  value  given  by  Moreno's  curves  [8]  to  two  digits 
(which  is  the  limit  of  accuracy  of  the  curves).  The  error  of  the  calculated 
value  from  the  measured  value  varies  from  0.9$  at  h/L=0  to  4.8$  at  h/L=5/6. 
The  true  error,  i .e.  between  the  simulation  and  the  lossless  cavity  which 
it  simulated,  was  probably  less  than  indicated  by  these  measured  values. 

The  reason  for  saying  this  is  that  the  reentrant  cavity  with  h/L=0  is  really 
a  simple  cavity,  and  the  theoretical  resonance  was  between  the  calculated 
and  measured  values.  The  calculated  value  was  only  in  error  by  0.2$  when 
compared  to  the  theoretical  value  at  h/L*0. 
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After  inserting  a  beam  hole  and  a  coupling  slot,  our  structure  is 
one  half  of  a  coupled  cavity,  whose  resonant  frequency  is  the  upper  cutoff 
of  our  structure.  (Actually,  this  is  one  half  of  an  "in-line"  slot  coupled 
cavity.  However,  the  change  in  resonant  frequency  of  a  single  cavity  in 
which  the  slots  were  in-line  as  compared  to  one  in  which  they  were  staggered 
was  measured  as  less  than  0 . 53£ . )  Figure  4.2  shows  the  measured  and  predicted 
(SOS)  resonance  values  at  various  stages  of  cavity  development.  It  is 
interesting  to  notice  that  the  accuracy  becomes  significantly  worse  when 
the  slot  is  introduced.  This  is,  in  fact,  the  point  at  which  the  problem 
truly  becomes  three  dimensional  (it  has  up  to  this  point  had  circular 
symmetry).  It  should  be  pointed  out  that  the  predicted  values  in  Figure  4.2 
are  actually  tiie  center  point  of  error  bars  (for  reasons  which  will  be 
explained  in  Chapter  V),  but  for  clarity  only  a  point  is  shown  at  this 
stage.  Although  the  predicted  upper  cutoff  value  is  in  error  by  approxi¬ 
mately  5%  (±  uncertainty),  we  are  getting  a  distinct  resonance  which 
encourages  us  to  try  to  find  more  points  on  the  Brillouin  curve. 

Simulations  of  multiple  cavities  were  performed,  which  predicted 
multiple  resonances  as  expected.  Two  cavities  with  the  grid  overlaid 
are  shown  in  Figure  4.3.  Notice  the  non-uniform  grid,  as  mentioned  in 
Chapter  II,  and  justified  in  Appendix  A.  These  multiple  cavity  runs  were 
all  done  with  a  ring  of  source  points  between  the  ferrules  in  one  cavity 
only,  oscillating  at  a  frequency  below  the  resonances.  Figure  4.4  is  the 
resulting  predicted  w-B.  The  solid  lines  are  drawn  through  the  physically 
measured  values.  If  the  best  end  of  the  error  bars  is  taken  as  the 
predicted  value,  then  the  worst  error  is  about  3?!.  Taking  the  worst  end 


of  the  error  bars,  the  predicted  value  Is  off  by  6%.  Therefore,  the  worst 
error  is  somewhere  between  3  and  6%.  Although  this  is  a  disappointingly 
large  value  for  the  error  in  prediction,  it  does  represent  a  slight  improve¬ 
ment  in  phase  velocity  prediction  over  typically  used  equivalent  circuit 
techniques.  In  Chapter  VI,  it  is  noted  that  by  application  of  a  "least 
squares"  cubic  to  the  data,  this  error  can  be  significantly  reduced.  For 
most  applications,  the  upper  and  lower  cutoff  are  not  the  areas  of  interest. 
The  area  in  which  we  would  desire  the  most  accuracy  would  be  centered  about 
Bp=tt/2,  say  from  n/3  to  2tt/3.  In  this  reduced  region,  using  the  least 
squares  curve,  the  error  is  everywhere  less  than  1$.  This  seems  to  be  a 
very  satisfying  result. 
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FIGURE  4.1  -  Measured  and  calculated  (by  the  code)  resonant 
frequencies  of  a  reentrant  cavity  with  various  values  of 
inner  conductor  height  (h).  All  dimensions  are  inches, 
the  two  radii  and  the  cavity  length  are  fixed,  the  inner 
conductor  height  being  the  only  variable. 
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FIGURE  4.3  -  Front  and  cross-section  views  of  a  two  cavity 
stack  with  the  grid  used  In  the  simulation  Indicated.  Notice 
that  the  grid  was  chosen  to  be  uniform  In  0  and  Z,  but  non- 
uniform  In  r. 


CHAPTER  V 
RESOLUTION 


Figure  5.1  shows  an  anatog  signal  which  is  just  resolved  by  the 
common  3  dB  definition.  If  discrete  samples  of  the  signal  are  taken,  as 
shown  in  Figure  5.1,  at  regular  intervals  of  fQ,  where  fQ  is  comparable 
in  magnitude  to  the  width  of  the  peaks,  then  we  see  (in  Figure  5.2)  that 
the  discrete  samples  are  not  resolved  by  the  3  dB  definition.  Furthermore, 
we  cannot  determine  from  the  discrete  samples  exactly  where  the  peak  is 
located.  We  can,  however,  say  that  the  peak  is  within  half  the  sampling 
width  of  the  local  maximum  sample.  With  these  two  factors  in  mind,  we 
have  defined  in  this  study  a  discrete  sample  to  be  resolved  if  there  exists 
a  sample  which  is  a  local  maximum.  That  is,  if  the  sample  is  greater  than 
both  of  its  nearest  neighbors.  Furthermore,  we  will  say  that  the  peak  is 
"resolved  to  fQ."  That  is. 


We  have  seen  that  the  spacing  in  frequency,  fQ,  is  determined  by 
the  total  run  time,  Nt0. 


Therefore,  the  total  run  time  should  be  long  enough  to  resolve  the 
resonances  in  two  ways.  First,  if  more  than  one  resonance  is  expected, 
then  fQ  must  be  small  enough  to  allow  both  of  the  peaks  to  be  local 
maxima.  Second,  fQ  must  be  small  enough  to  determine  the  location  of 


each  peak  to  reasonable  accuracy.  One  would  expect  then,  that  the  limiting 
factor  in  resolution  would  be  the  cost  associated  with  run  time. 

The  limiting  factor  for  this  study,  turned  out  to  be  an  instability 
in  the  calculations.  For  an  undetermined  reason,  it  was  found  that  the 
higher  frequencies  grew  in  amplitude  much  faster  than  the  low  ones.  When 
these  high  frequencies  became  extremely  large  in  amplitude,  the  lower 
frequencies  (which  were  in  the  range  of  interest)  were  obscured  by  them. 
Figures  5.3  through  5.5  illustrate  this  problem.  In  Figure  5.3,  the  source 
at  3.0  GHz,  and  the  lowest  resonance,  expected  just  below  4.0  GHz,  are  not 
yet  resolved  due  to  the  large  size  of  the  frequency  spacing  (0.4  GHz). 

After  doubling  the  run  time,  Figure  5.4  shows  that  the  source  and  primary 
resonance  are  resolved,  but  most  of  the  energy  is  now  in  the  higher 
frequencies,  with  the  major  resonance  at  about  46  GHz.  Figure  5.5  shows 
that  after  again  doubling  the  run  time,  the  source  and  primary  resonance 
are  now  totally  obscured  by  the  resonance  at  46  GHz. 

Since  this  instability  was  one  which  grew  in  time,  it  was  possible 
to  perform  one  long  run  and  take  the  data  for  the  FFT  from  time  zero  up 
until  the  instability  just  obscured  our  desired  results  and  then  truncate 
all  following  data.  This  was  the  method  used  to  obtain  the  w-B  plot  shown 
in  Chapter  V.  In  order  to  obtain  the  data  shown  there  for  the  160°  slot, 
two  runs  were  performed  which  modeled  three  and  four  cavities,  and  gathered 
slightly  more  data  than  could  be  used  (due  to  the  instability).  These 
runs  were  performed  at  a  total  cost  of  less  than  $150. 


ANALOG  SIGNAL  WITH 
DISCRETE  SAMPLINGS 


FIGURE  5.1  -  Discrete  samples  overlaid  on  an  analog  signal 
with  two  peaks  that  are  just  resolved  by  the  common  3  dB 
def in  it  ion. 
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DISCRETE  REPRESENTATION 
OF  THE  SIGNAL 


FIGURE  5.2  -  Discrete  samples  from  Figure  5.1.  The  samples  are 
not  resolved  by  the  3  dB  definition,  and  furthermore,  the  true 
location  in  frequency  of  the  maximum  may  only  be  determined  to 
the  limits  indicated. 
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FIGURE  5.4  By  doubting  the  sampling  rate  (to  every  200  MHz) 
the  resonance  and  source  are  resolved,  but  the  high  frequency 
amplitudes  have  grown. 
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FIGURE  5.5  -  After  doubling  the  sampling  rate  again  (to  every 
100  MHz)  the  enormous  amplitude  at  high  frequency  has  obscured 
the  area  of  interest. 
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CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 


1+  has  been  pointed  out  to  the  author,  that  the  shape  of  the  predicted 
co-8  is  quite  unlike  the  curves  which  are  characteristic  of  staggered  slot 
coupled  cavity  circuits  (first  derivative  zero  at  Bp=0  and  tt,  and  the 
second  derivative  changing  signs  at  approximately  8p=tt/2).  While  this 
objection  is  true,  it  is  a  fact  that  may  be  used  to  one's  advantage.  By 
applying  a  least-squares  fit  to  a  cubic  function  with  zero  first  deriva¬ 
tives  at  the  end  points  (see  Appendix  D),  it  is  possible  to  determine  a 
curve  which  is  fairly  close  to  the  measured  one.  Figure  6.1  demonstrates 
this  method,  where  the  center  points  of  the  error  bars  has  been  used  for 
the  data  in  equation  C. 2.  The  predicted  curves  look  very  much  as  expected 
now.  The  maximum  error,  which  occurs  at  the  upper  cutoff  of  the  160°  slot, 
is  now  about  3.9$. 

While  this  project  has  suffered  from  some  problems,  especially  in 
getting  a  desirable  level  of  resolution,  we  have  obtained  our  primary 
objective.  That  is,  we  have  shown  that  an  existing  plasma  simulation 
code  may  be  used  to  predict,  to  a  fair  degree  of  accuracy,  the  cold  cir¬ 
cuit  characteristics  of  a  microwave  circuit. 

A  very  promising  result  is  that  the  accuracy  was  not  cost-limited, 
but  rather  was  limited  by  the  stability  problem  as  explained  in  Chapter  V. 
This  is  promising  because  if  we  were  cost-limited,  we  would  be  at  the 


mercy  of  the  computer  industry.  That  is,  no  further  significant  progress 
could  be  made  until  a  faster,  more  efficient  computer  was  built.  As  is, 
there  is  good  hope  of  solving  the  problem  and  going  on  to  improved 
solut ions. 

We  give  now  some  cost  figures  to  illustrate  the  typical  expense  of 
this  method.  The  two  runs  (three  and  four  cavities)  needed  to  find  the 
data  for  the  160°  slot  shown  in  Figure  6.1  were  performed  at  a  combined 
cost  of  under  $150.  These  runs  were  performed  overnight  in  batch  mode 
on  a  slow  que  (low  priority),  which  results  in  a  charge  of  $12  per  CPU 
hour.  Therefore,  the  two  runs  required  about  12  CPU  hours  on  the  VAX 
11/780.  This  machine  has  a  speed  of  approximately  0.5  Mflop. 

Future  investigation  should  yield  solutions  to  the  high  frequency 
instability.  The  solution  may  be  a  combination  of  changing  source 
methods,  grid  location,  and  applying  low-pass  filters  to  the  fields  in 
the  numerical  algorithms. 

Subsequent  to  the  completion  of  this  project,  the  high  frequency 
noise  problem  experienced  by  the  author  was  verified  in  a  run  performed 
by  MRC  personnel.  It  should  be  noted  that  this  problem  is  unique  to 
the  cylindrical  coordinate  version  of  the  code.  The  rectangular  version 
has  been  in  existence  for  some  time,  and  due  to  numerous  test  cases,  there 
is  a  high  level  of  confidence  in  its  accuracy;  Further  work  is  currently 
being  done  at  MRC  to  determine  the  source  of  this  problem  and  to  alleviate  it 
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FIGURE  6.1  -  Measured  and  calculated  values  for  the  Brillouin 
plot  of  the  coupled  cavity  circuits.  By  applying  a  least  squares 
fit  to  the  calculated  data,  we  see  that  it  agrees  quite  well  with 
the  measured  data.  In  fact,  over  the  middle  one  third  of  the 
band,  which  is  the  area  of  interest,  the  agreement  is  everywhere 
within  1i  error. 
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APPENDIX  A 


A  DERIVATION  OF  THE  SECOND-ORDER  ACCURACY 
OF  THE  CENTERED-DIFFERENCE  DERIVATIVE 
FOR  A  NON-UNIFORM  GRID 


This  appendix  shows  that,  under  certain  conditions,  the  centered 
difference  approximation  to  the  derivative  of  a  function  will  be  accurate 
to  second  order  without  the  restriction  of  a  uniform  grid.  The  signifi¬ 
cance  of  this  fact  is  that  it  demonstrates  both  the  flexibility  and  the 
integrity  of  the  field  algorithms  used  in  this  study. 

The  centered  difference  formula  is  given  by 

f(Xj*l  )  — f ( x | _i ) 

F'Cxs)  =  - ' - 

where  f  denotes  the  true  function  (whose  values  are  known)  and  F'  denotes 
the  approximation  to  the  true  derivative  f’  (whose  value  is  not  known). 
Now,  if  all  derivatives  of  f  exist  on  the  interval  Xj_1<x<xj+l>  then  f  may 
be  expanded  in  a  Taylor  series  by  defining 

h  j  =  x  j -x  j  _j 

so  that 

xi-  i**rhi 


xi  +i=x  i+*1i  +  i  • 


Then  the  expansion  is 


h-2  h-3 

f  (  x  j  — h  j  )  =  f  ( x  j  )-hjf’(xi  )*  -j-  f"(xi  J-yp  f’"  <xj  >♦... 


tv2  h1  * 

f(xj+h;+j )  -  f  (  x  ,  )  +h  j  *  j  f  '  (  x  j  )  ♦  f"(x|)  ♦  -ilL.  f  ” » ( x  I )  ♦ . 


and  the  numerical  derivative  is 
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F  *  (  x  i ) 


2  .  2, 


=  hy--h7  +  )f  *(xi  )*’>i(hWl-hi  JfCxj ) 
+  37  (hj3+i+hi3)f"'  (*i >  +  .  •  •> 

=  f  * ( x  j ) *\ ( h  j  +  j -h  j )  f " ( x  j )  ■ 


(hu-hj-}  f,M  ( x ;  . 


6(hi  +  1+hn 


(A.  1 ) 


From  which  it  is  obvious  that  second  order  accuracy  [0(h2)]  is  obtained 
for  the  uniform  spacing  case;  i.e.  hj+1=  hj.  However,  if  we  require  the 
Xj  to  be  a  function  of  a  uniformly  spaced  parameter  n, 

Hj  =  r|0+(i“1)^T1  =  constant  (A. 2) 

then  second  order  accuracy  may  be  retained  without  the  stringent  requirement 
of  a  uni  form  grid. 

Proof: 

The  functional  dependence  may  be  expressed  as 

xrx0  =  a  ( n  j  -nQ )  ( n  j  _n0  > 2  (A,3) 

and  ( A . 1  )  may  be  rewritten  as 

F  *  ( x  j  )  =  f  *  ( x  j  )+vj(  x  j  +  1-2x  j+xi_1  )f"(xj  )  +  . . . 

=  f  ’  ( X i  )-*M 0<o+a(ni  +  ,-r)0)+b(nj  +1-n0)2J 
-2[x0+a(ni-r)o)+b(rii-no)2]  +  [xo+a(T1i-i"rlo)+b(ni-rT10)2]}f"(xi )  +  -*- 


-41- 


?•] 


9 

I 


Now  we  make  use  of  (A. 2) 

n  j  + 1  _r)0  =  i6n;  n  j  -n0  =  c  1-1  )6n;  n ;  _  x  -n0  =  n-2)6n 

3o  that  the  numerical  derivative  reduces  to 

F'(xj)  =  f'(x,  )+H{[a(  i6n)+b(  i6n)2]-2[a(  i-1  )6n+b(  i-1  )2(6n)2J 
+  (a ( i-2)6n+b( i-2)2(6n)2]}f"(xi )  +  ... 

=  f '  ( x  j  )+wj{a6n[i-2(  i-1 }  +  ( i-1 )  +  ( i-2)]+b6n2[i2-2(  i-1  )2  +  (  i-2)2]} 


f  * ( x ] )+b6n2f"(x j )  +  . . 


T  hus , 


F '  ( x ;  )  =  f * ( x ; )+0(6n2) 


and  the  numerical  derivative  is  accurate  to  second  order. 
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APPENDIX  B 

THE  DISCRETE  FOURIER  TRANSFORM 


This  appendix  formulates  the  DFT  pair  which  is  used  in  this  report 
to  gain  frequency  information  from  a  limited  time  history  of  the  fields. 
Although  these  equations  may  be  found  in  texts  on  signal  processing  [5], 
it  seemed  well  to  reproduce  them  here  to  clarify  the  relation  between  the 
generalized  parameters  and  the  parameters  of  interest  to  us,  which  are 
time  and  frequency.  This  also  shows  the  tradeoff  between  computer  cost 
(total  number  of  samples)  and  resolution  (spacing  in  frequency  samples). 

In  general,  a  function,  f(t),  on  the  interval  0  <  t  <  T  can  be  expressed 
as  an  infinite  sum  of  exponentials 


f(t)  =  Z  Cn  ej(2n*/T>t 

n  =  -® 


If  the  function  is  sampled  at  discrete  points  such  that 


t  El 

N 


m  =  1 ,  2,  3,  . . . ,  N 


then 


f(t)  =  f(m)  =  Z  cneJ(2n1T/N)m- 


n  =  -°o 


Now,  we  can  take  advantage  of  the  periodicity  of  the  exponential 
by  substituting  K'^N  for  n,  where  K  takes  on  values  only  in  the  range  from 
1  to  N.  Then  we  have 


f(m)  =  Z  Z  CK  £  e j( 2Ktt/N )m  ej( 28-Ntt/N )m 
£=-®K=1 


=  Z  Z  C«,Jl  eJ(2KlT/N)rn 
Z—co  K=  1 


N  00 

=  Z  ej(2KTT/N)m  Z  CK  l 

K= 1 
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The  last  term  is  merely  a  summation  of  constants,  so  we  can  define  a 


new  constant 


or  equivalently 


F(K)  =  2  Ck.jl  . 

£.  =  -» 


f (m)  =  2  F(K )  2l<lt^ )m 

K  =  1 


f (m )  =  2  F(n)ej(2mr/N)m  . 

n=0 


In  order  to  evaluate  the  coefficients  F(n),  we  multiply  both  sides 
by  e~J2Tr®'rr,/^  and  sum  from  m=0  to  N-1. 

N-1  N-1 

2  f(m)e-j2iTilm/N  =  2  £  F(n)ej2iTm(n-fc)/N  . 

m=0  m=0  n=0 

At  this  point,  we  make  use  of  the  orthogonality  relation 


(B.1) 


N-1 

l  ej2TTmr/N 
m=0 


1  r=kN 
0  otherwise 


k  an  integer 


By  letting  n-£=r  in  (B.1),  we  see  that  the  RH5  is  only  non-zero  if 

n=£  +  r=5.+kN.  Since  n=0,  1,  _ ,  N-1  this  is  only  true  at  one  point.  For 

convenience,  let  8,  be  in  the  range  0  to  N-1,  then  k  must  be  zero  for  the 
condition  to  be  satisfied.  Thus, 

N-1  N-1 

2  F ( n )  2  eJ’2  m(n-H)/N  =  Fa).N 

n=0  m=0 
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APPENDIX  C 


A  LEAST  SQUARES  FIT  TO  A 
CUBIC  POLYNOMIAL 


The  general  equations  for  a  least  squares  fit  of  a  discrete  set  of 
data  points  to  a  polynomial  of  any  order  may  be  found  in  most  texts  on 
numerical  analysis  [9].  However,  with  the  inclusion  of  boundary  conditions 
on  the  curve,  some  of  the  variables  may  be  eliminated  and  the  equations 
are  changed.  For  our  problem,  the  value  of  the  curve  at  the  endpoints 
is  not  known,  but  the  first  derivative  is.  This  appendix,  therefore, 
develops  the  equations  necessary  to  find  the  cubic  curve  which  minimizes 
the  square  of  the  error  from  the  data  points  while  retaining  a  zero  first 
derivative  at  the  endpoints. 

Given  a  set  of  data  points  [xi;yf]  1=1,  2,  . . . ,  m  on  the  interval 
0<x<l  (note  that  all  intervals  may  be  normalized  to  this)  and  the  cubic 
i nferpolanf 

Y=a+bx+cx2+dx3 

dY 

and  the  constraints  -y~  =  0  at  the  endpoints.  We  wish  to  minimize  the  error 


m 

E  =  I  (Yj-y ; )2 
i  =  1 

with  respect  to  the  variables  a,  b,  c,  and  d.  However, 

two  of  the  variables  because  of  the  constraints. 

Y'  =  ~  =  b+2c+3dx2 
dx 

Y'(O)  =  b  =  0 

and 


So  that 


Y'  1)  =  2c-*  3d  =  0 


Y  =  a*dx2  (x~4) 


we  can  el imi nate 
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Now  we  minimize  the  error  with  respect  to  the  two  variables,  a  and  d. 


||  =  2  (Yj-yj)  =  0 


i  =  1 


dd  i=1  3d 


but  we  see  from  (C.1)  that 
3Yj 


=  1 


3  a 


3  Y ;  «  x  j 2 ( x  j  —3/2 ) 
3d 


therefore  our  equations  are 
m 

Z  (Yj-y j )  =  0 
i  =  1 


Z  ( Y  j -y i ) x  j 2 ( x ; —3/2 ) 
i  =  1 


V  ( C .2  ) 


=  0 


.y 


where 


Yj  =  a+dxj2(xj-3/2). 


[  So  we  have  two  equations  from  which  our  two  unknowns  (a  and  d)  may  be  found. 
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